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At large Nc, cold nuclear matter is expected to form a crystal and thus spontaneously break 
translational symmetry. The description of chiral symmetry breaking and translational symmetry 
breaking can become intertwined. Here, the focus is on aspects of chiral symmetry breaking and its 
possible restoration that are by construction independent of the nature of translational symmetry 
breaking — namely spatial averages of chiral order parameters. A system will be considered to be 
chirally restored provided all spatially-averaged chiral order parameters are zero. A critical question 
is whether chiral restoration in this sense is possible for phases in which chiral order parameters 
are locally non-zero but whose spatial averages all vanish. We show that this is not possible unless 
all chirally-invariant observables are spatially uniform. This result is first derived for Skyrme-type 
models, which are based on a non-linear sigma model and by construction break chiral symmetry 
on a point-by-point basis. A no-go theorem for chiral restoration (in the average sense) for all 
models of this type is obtained by showing that in these models there exist chirally symmetric order 
parameters which cannot be spatially uniform. Next we show that the no-go theorem applies to large 
Nc QCD in any phase which has a non-zero but spatially varying chiral condensate. The theorem 
is demonstrated by showing that in a putative chirally restored phase, the field configuration can 
be reduced to that of a nonlinear sigma model. It is also shown that this no-go theorem is fully 
consistent with the vanishing of the spatial average of the chiral condensate ^Tr(f/) (as happens 
in "half-skyrmion" configurations). This is because the chiral condensate is only one of an infinite 
set of chiral order parameters, some of which must be non-zero. It is also shown that while an 
approximation of a unit cell of a Skyrme crystal as a hypersphere does lead to a phase which is 
chirally restored (in the average sense), this is an artifact of the approximation. 

PACS numbers: 11.15.Pg, 12. 38. Aw, 21.65.-1 



I. INTRODUCTION 

One of the long-standing problems in QCD is to un- 
derstand finite density nuclear matter. Such matter is 
only poorly understood in terms of QCD but is of fun- 
damental importance in developing a theoretical under- 
standing of the qualitative and quantitative features of 
QCD equations of state and phase diagrams. It is also 
important in studying astrophysical objects such as neu- 
tron stars, which are composed of dense nuclear matter, 
and in heavy ion collisions. While zero density nuclear 
matter can be studied using lattice QCD, at finite den- 
sities or equivalently chemical potential and low temper- 
atures, the fermion sign problem renders lattice studies 
intractable. 

Given this situation one is often forced to rely on Q CD- 
inspired models to get insight. This is problematic. First 
of all, it is obvious that the models are not QCD and one 



'Electronic address: 'prabal@umd.eduj 

t Electronic address: cohen@physics . umd.edu| 
■f Electronic address: imcstro@umd.edul 



needs to be cautious in interpreting their results. Sec- 
ondly, many models rely at least implicitly on QCD be- 
ing close to the large Nc limit [l|, 0| where mean-field 
methods can be justified. At large N^, quantum effects 
become negligible and nuclear matter behaves classically 
in some important sense. In this limit of high densities 
and large Nc, it is generally believed that it is energet- 
ically favorable for nuclear matter to crystallize in its 
ground state. 

The large Nc limit was recently used to motivate the 
existence of a new phase of finite density nuclear mat- 
ter. McLerran and Pisarski Q argued that at low tem- 
peratures nuclear matter condenses into a "quarkyonic 
phase" of confined but chirally-restored nuclear matter 
dominated by baryon-baryon interactions. This idea was 
further developed in Ref. Ql- This phase is supposed 
to occur in large Nc QCD in the regime in which the 
quark chemical potential is of order iV^ large compared to 
Aqcd • The logic is simply that in this regime the quarks 
do not effect the gluons and thus confinement described 
in terms of the Polyakov line is unaffected. Thus, pro- 
vided the chemical potential is large enough to cause chi- 
ral restoration, one necessarily enters a quarkyonic phase. 
This analysis is, strictly speaking, only clean in the chiral 



2 



limit of zero quark masses and for the purposes of this 
paper, we will ignore the effects of quark masses and set 
them to zero at the outset. 

There is a potential difficulty with the argument Q. 
A priori it is not trivially obvious from the structure of 
large Nc QCD that at high baryon density, the system 
necessarily becomes chirally-restored. This issue is at 
the crux of this paper. The primary question is whether 
chiral symmetry is restored at high density. We have 
no answer to this critical question. However, there is a 
secondary question of importance which we can address, 
namely if chiral symmetry is restored, by what mecha- 
nism does this occur? 

To set the context for this, let us consider a model- 
dependent argument used in Ref. Q to justify the as- 
sumption that chiral symmetry is restored in large Nc 
QCD at sufficiently high density. The argument is based 
on the Skyrme model |^ treated at the classical level 
(which is justified at large Nc). The model, although 
phenomenological, does capture many aspects of large 
Nc QCD and of chiral symmetry, ft has proven to be 
a reasonable, if crude way to describe static properties 
of baryons such as masses, charge radii, and magnetic 
moments Q. It has been used to describe nuclear mat- 
ter which at mean-field level is naturally described as a 
crystal 

At first sight, one might think that the Skryme model, 
which is based on a non-linear sigma model and hence 
by construction breaks chiral symmetry at every point 
in space, is unsuitable to study chiral restoration. How- 
ever, at this stage, it is worth noting that at low but 
non-zero baryon density, two distinct types of symmetry 
are spontaneously broken: chiral symmetry, which is bro- 
ken as in the vacuum, and translational symmetry, which 
is broken by the formation of a crystal. The two types of 
spontaneous symmetry breaking can become intertwined. 
Since our focus, is on chiral symmetry breaking and not 
on translational symmetry breaking it is reasonable to 
consider only those aspects of chiral symmetry breaking 
which do not depend on translational symmetry break- 
ing. 

It is easy to construct order parameters that are sen- 
sitive to chiral symmetry breaking but insensitive to any 
details of translational symmetry — namely by spatially 
averaging standard chiral order parameters. It is not un- 
reasonable to consider a system chirally-restored if all 
of the spatially averaged chiral order parameters vanish; 
this paper will focus on the question of chiral restoration 
in this spatially averaged sense. In the remainder of the 
paper, the phrase "chiral restoration" will often be used as 
a shorthand for "chiral restoration in the spatial averaged 
sense". One might hope that in Skyrme-type models, chi- 
ral symmetry, although broken at every point in space, 
could be restored in this spatially averaged sense. 

In arguing that large Nc QCD should restore chiral 
symmetry at sufficiently high density, Ref. Q presents 
evidence that the Skyrme model treated classically, in 
fact, exhibits chiral restoration in this spatially-averaged 



sense and argues that large Nc QCD might similarly ex- 
hibit chiral restoration since the Skyrme model captures 
many aspects of large Nc QCD. The evidence that the 
Skyrme model restores chiral symmetry at high density is 
of two types: (i) the vanishing of the spatially-averaged 
chiral condensate found in studies of different types of 
Skyrme crystals (7l-[l0l. [l^ . [isl . [Tsl and (ii) the vanishing 
of spatially-averaged chiral order parameters as seen in 
analytical studies in which a Skymion in the crystal is 
approximated as a single Skyrmion on the compactified 
manifold of a hypersphere [8j, ll6Hl8| . 

The evidence that chiral symmetry is restored in the 
spatially-averaged sense at sufficiently high density in the 
Skyrme model appears to be quite strong. However, as 
will be shown in this paper, it is misleading. We prove a 
theorem that it is impossible for any Skyrme-type model 
to restore chiral symmetry in the spatially-averaged sense 
at non-zero baryon density. We further show that the 
no-go theorem also applies to large Nc QCD with a non- 
zero, spatially varying chiral condensate. This is done 
by showing that in a chirally-restored phase, expectation 
value of operators can be reduced to that of a generic 
nonlinear sigma model. 

Before proving this no-go theorem, it is useful to dis- 
cuss how it can possibly be correct given the evidence 
presented in Ref. @. First, let us consider the evidence 
of the vanishing of spatially-averaged chiral condensate 
in Skyrme crystals above some critical density. It is cer- 
tainly possible that the chiral condensate itself may van- 
ish on spatial integration and indeed, it is expected to 
happen in so-called "half-skyrmion" configurations of the 
crystal's unit cell (sl. [Tll - flSl . [l5| . This occurs when iso- 
lated skyrmions (whose interaction energy at low den- 
sity scales as r^) [T3|) are brought close together. It is 
believed that lowest possible energy consistent with the 
topology of the model (a winding number of unity per 
unit cell of the crystal) is a new type of phase based on 
a configuration that has an additional axial symmetry. 
This constrains the field to have a — —1 at the center 
of the unit cell. The configuration has a baryon number 
of i and can be transformed into another half-skyrmion 
by the transformation (cr, tt) (— cr, — tt). The unit cell 
consists of two of these half-skyrmions. The chiral order 
parameter ^Tr{U) = a, where cr (which is assumed to 
be proportional to qq of QCD), vanishes when integrated 
over the unit cell. 

The key point is that while chiral restoration (in the 
average sense) necessarily implies the vanishing of the 
spatially-averaged chiral condensate, the converse is not 
true. It is logically possible that chiral condensate could 
vanish due to the restoration of a discrete symmetry while 
other spatially-averaged chiral order parameters remain 
nonzero. A priori it may seem implausible that such 
a scenario could be realized. Indeed, it is reminiscent 
of Stern's suggestion that in the vacuum, (qq) vanishes 
but chiral symmetry remains broken due to higher di- 
mensional condensates [l9j . However, as shown by Ko- 
gan, Kovner and Shifman (KKS) [20| using an elegant 
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argument based on a rigorous QCD inequality for the 
Euclidean space functional integral, the Stern scenario 
is inconsistent with QCD. One might worry whether a 
similar scenario is also impossible for finite density mat- 
ter. However, the KKS argument depends on a real and 
positive functional determinant, which while true for the 
vacuum does not hold for finite density matter in which 
the chemical potential leads to a functional determinant 
that is not manifestly real and positive. Thus, this argu- 
ment is not applicable to finite density matter. Moreover, 
as will be shown in this paper, this unlikely seeming sce- 
nario must occur for both Skyrme models and for large 
Nc QCD if there exists a regime where the average chiral 
condensate is not zero locally but vanishes under spatial 
averaging. 

Now let us consider the second class of evidence for 
chiral restoration in the average sense — namely the fact 
that when a Skyrmion in the crystal is approximated as 
a single Skyrmion on the compactifled manifold of a hy- 
persphere, all spatially-averaged chiral order parameters 
vanish. In effect by using Skyrmions placed on the com- 
pactifled manifold of a hypersphere "surface", S^{L) em- 
bedded in [l^ [13, HH , one is effectively chopping 
off the "edges" of a unit cell in a real Skyrme crystal 
and somehow mapping this onto the hypersphere again 
with the winding number set to one. Of course, in do- 
ing this one is making an uncontrolled approximation, 
but the virtue is simplicity with many quantities ana- 
lytically calculable. The density of the Skyrme crystal 
is assumed to be the inverse volume of the hypersphere. 
For low densities (or large radius hypersphere) , there are 
two solutions — a stable low energy one in which the wind- 
ing number density (i.e., baryon density) is localized and 
spatially-averaged chiral order parameters are nonzero 
and an energetically unstable phase in which the winding 
number density is spread uniformly and spatial averages 
of all chiral order parameters vanish. However at some 
critical density (or equivalently radius) the two solutions 
merge and above this density the lowest energy solution 
is the one with uniform baryon density and chiral restora- 
tion in the spatially-averaged sense. 

Of course, this evidence for chiral restoration depends 
on an uncontrolled approximation. One does not know 
at the outset how well one can approximate a Skyrmion 
in the crystal by a single Skyrmion on a hypersphere. 
Presumably, the logic underlying the approximation is 
that the principal effect of putting a Skyrmion into a 
crystal is to restrict the space over which it can spread. 
If this is the case it is natural to assume that a Skyrmion 
in the restricted space of a hypersphere is likely to be 
qualitatively similar to a Skyrmion in a crystal. Thus, it 
is plausible that while details of such a precise equation of 
state or the density of the phase transition to the chirally- 
restored phase would be affected by the approximation, 
qualitative features such as the existence of a chirally- 
restored phase would not. 

The preceding argument depends on a critical assump- 
tion, namely that using a hypersphere to restrict the vol- 



ume of the Skyrmion acts generically like other restric- 
tions on its volume. This need not be true. The hy- 
persphere is a rather special geometrical structure. One 
might worry that the chiral restoration seen on the hy- 
persphere is not a generic feature at all, but rather an 
accidental feature of the ad hoc choice of geometry. If 
this is the case then the fact that Skyrmions exhibit chi- 
ral restoration at high density on the hypersphere gives 
no insight into what happens for Skyrme crystals. As 
will be shown, the vanishing of all spatially-averaged chi- 
ral order parameters is a result of the peculiarities of the 
geometry of the hypersphere. This can be seen clearly 
by keeping the Skyrmion on a compact space but one 
in which the geometry is distorted away from the hy- 
persphere. When this is done, chiral restoration (in the 
average sense) disappears. 

The paper is organized as follows. In the following sec- 
tion, we prove a no-go theorem that prevents chiral sym- 
metry restoration (in the average sense) for all Skyrme- 
type models, i.e., non-linear sigma models in which the 
winding number density is identified as the baryon den- 
sity. This proof consists of two parts. First, we derive 
as a necessary condition for chiral restoration (in the av- 
erage sense) the condition that all chiral singlet observ- 
ables must be spatially uniform. Next, we show that 
a particular chirally symmetric order parameter cannot 
be uniform while chiral symmetry is restored simulta- 
neously in Skyrme-type models, thus establishing the 
theorem. Then we extend our no-go theorem from the 
realm of models to large Nc QCD itself. In particular, 
we prove that large Nc QCD at non-zero baryon density 
cannot have a chiral condensate, which is non-zero but 
non-uniform spatially in such a way that all spatially- 
averaged chiral order parameters vanish. In the follow- 
ing section, we discuss in detail the case of a Skyrmion 
confined to a hypersphere in light of the no-go theorem. 
We show explicitly why the proof of the no-go theorem 
works for crystals in flat space but breaks down for a 
Skyrmion confined to a hypersphere. Next, we discuss 
why the insights gleaned from a Skyrmion confined from 
a hypersphere is misleading. We do this by consider- 
ing Skyrmions confined to a compact manifold which is 
deformed from a hypersphere; chiral symmetry in the av- 
erage sense is not restored. This demonstrates that chiral 
restoration observed on the hypersphere is an artifact of 
the (unphysical) choice of geometry. We end with a brief 
discussion of these results. 



II. A NO-GO THEOREM IN THE SKYRME 
MODEL 

In this section, we prove a no-go theorem for Skyrme- 
type models, which we define as chirally invariant non- 
linear sigma models based on a matrix- valued field: U (x) 
with U £ SU{2) (and possibly other fields) with a La- 
grangian rich enough to support stable topological soli- 
tons. The model is treated classically, in keeping with 
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large Nc QCD. The classical field, U{x), which we take 
to be continuous is a mapping from — >■ SU{2) for flat 
space (and from S^{L) — > SU{2) for skyrmions on the hy- 
persphere). Physically, the field U encodes the dynamics 
of the pion. Such models automatically have an alge- 
braically conserved winding number current, (that is 
dfj^wf^ ~ for any field configuration) where 

,A.-/? Tr [(c/ta,[/)(c/ta„c/)(c/ta^[/)] 

^ = 2A^2 ("■!) 

In these models, is often taken to be to equal to the 
conserved baryon current, although in principle they can 
differ by a quantity, which is conserved and has zero net 
charge. 

In proving the no-go theorem, the precise details of the 
Lagrangian play no role nor does the identification of the 
winding number current with baryon current. However, 
its chiral properties are essential. The Lagrangian must 
be invariant under a global SU {2)l x SU{2) h transforma- 
tion where [/(x) ^ LU{^)R) and C/t(x) ^ RU\^)L'^ . 
If the theory has other fields, they must transform in a 
consistent way. Of course, the non-linear sigma mod- 
els build in spontaneous symmetry breaking at the out- 
set; while the Lagrangian is chirally invariant, any field 
configuration at any point in space is not. That is 
C/(x) 7^ i?[/(x)it. Thus, by construction, chiral sym- 
metry when evaluated point-by-point in space cannot he 
restored in any model in this class. 

The notion of chiral symmetry restoration in the 
Skyrme model only makes sense as a global property and 
not a local one. We define chiral symmetry to be restored 
in the spatially-averaged sense if all functions of the clas- 
sical fields when spatially averaged over a unit cell equals 
the average of the same function averaged uniformly over 
the internal space of chiral rotations. One particularly in- 
teresting class of observables are those that depend on U 
(which already encodes spontaneous symmetry breaking) 
but not on spatial derivatives (or other possible fields in 
the problem) where one can project out the chiral singlet 
part by averaging uniformly over the SU{2) field. Thus, 
for example, chiral restoration in this sense requires that 
any scalar function of U (x) only (and not its derivatives) 
will satisfy 

f dVF{U{^))^^[ d^iF{U{ti)), [11.2) 
Ji?3 27r Jsu{2) 

where F is an arbitrary real-valued function of U and 
/i represents the three angles needed to specify a general 
SU{2) matrix and d^ is the Haar measure of SU{2). The 
logic of this is quite simple: averaging over the internal 
space ensures that all chiral order parameters will vanish 
while chiral sing lets wiU not. Note that in Eq. (|II.2|1 the 
average of the function of the field U^x) in the internal 
space is independent of the dynamics of the model, but 
its spatial average depends on the equations of motion. 

In this section, we prove that chiral restoration in 
the spatially-averaged sense is impossible for the Skyrme 



model. This proof has two parts. The first is a demon- 
stration that chiral restoration in this sense cannot oc- 
cur unless all chiral singlet observables are spatially uni- 
form. The second is that chiral restoration is inconsistent 
with the fact that all chiral order parameters must van- 
ish (upon spatially averaging) . This will be demonstrated 
via the explicit construction of a chiral singlet operator 
which cannot be spatially uniform in a regime where the 
spatial average of all chiral order parameters vanish. 



A. Chiral restoration in the average sense requires 
spatial uniformity for all chiral singlets 

The key to the first part of the proof is the intuition 
that the only natural way chiral restoration in the average 
sense can occur is if spatially integrating corresponds to 
integrating over the internal space with uniform weight- 
ing. The central part of the proof is a demonstration 
that this intuition is correct. Ultimately, this proves to 
be such a strong constraint that no field distributions in 
Skyrme-type models can satisfy it unless all chiral singlet 
observables are spatially uniform. 

To begin the formal treatment of the first part of the 
proof, we note that chiral restoration requires all chiral 
order parameters vanish, or equivalently that the spatial 
average of all functions equals the uniform average over 
all chiral rotations. Thus, to show that chiral restora- 
tion cannot occur it is sufficient to show that there exists 
some subset of chiral order parameters for which it is 
not true. Here, we will focus on observables constructed 
entirely from the field [/, which is local and hence de- 
pends on a single point x, as in Eq. (jll.2[) . Examples 
of this class of observables include F{U) = Tr[[/], and 
F{U) = Tr[J7Ti]2 -hTr[;7]2. Note that function U is com- 
pletely specified by three Euler angles, denoted collec- 
tively as ^ so that the most general F is simply the most 
general function ^ consistent with periodicity conditions. 
The Wigner D-matrices provide a complete basis of func- 
tions, which satisfy these boundary conditions so that we 
can always decompose F {U{ijl)) into the following form: 

where the coefficients f^^ ^, are constants which depend 
only on the function F . Note that the mapping from 
physical space onto U can be recast as a mapping from 
real space onto /^: F (J7(x)) = F (C/(/^(x))) We will show 
that for any smooth field configuration, there must be 
members of this class of observable which do not satisfy 
Eq. (|II.2[) unless all chiral singlet observables are spatially 
uniform. 

To proceed, it is useful to recall that the Wigner D- 
matrices satisfy an orthogonality condition 

/ dnDi,^^,D*l^ = ^^3^,^6r,5^'k . (II.4) 
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Since, the D-matrix for j = is a constant independent 
of /i, it follows from the orthogonality relation that for 



(II.5) 



Thus, it is clear that is a chiral order parameter (i.e., 
an observable which is necessarily zero if chiral symmetry 
is unbroken) if, and only if, when written in the form of 
the decomposition of Eq. pi.3[) the coefficient /q o = 0- 

Since the Wigner-D matrices form a basis, Eq. (jll.2[) 
can only be satisfied for all choices of F if it is satisfied for 
each Wigner-D matrix. A simple way to implement this 
is to simply choose F = j, (j ^ 0). Equation (|II.4p 
then implies that if chiral restoration in the average sense 
occurs, then 



F{u{^l{^)))dv 



Dl.kJ{l^)dl^^^. (11.' 



where J{^) is the determinant of the Jacobian for a trans- 
formation from the physical space onto the internal space 
and the right-hand side of Eq. (jll.4[) is zero by orthogo- 
nality provided that j 7^ 0. Note that the mapping of x 
onto /i need not be one-to-one but one can always cast 
the integral into the preceding form by combining the Ja- 
cobian determinants for any regions in real space which 
are mapped into the same region of internal space. 

The next step is simply to note that since the Jaco- 
bian determinant is a scalar function of it too can be 
decomposed into a basis of Wigner-D matrices 



Ji.1^)^ E 



C^m'.k'Kl'k'il^)- 



(II.7) 



j' ,7n' .k' 



Here, we have used complex conjugates of the Wigner 
matrices to exploit the orthogonality condition of Eq. 
()II.4[) . Note that the coefficients c;^, characterize the 
field configuration and are independent of the observ- 
able F. Inserting this decomposition into Eq. (jII.6P for 
F ~ D^^ ^ and using orthogonality implies that if chiral 

restoration in the average sense occurs only if c^„ = 
for all J 7^ 0. 

Note that the choice of F is arbitrary and thus we we 
can repeat the analysis for all j, m and k and thus de- 
duce that if chiral symmetry is restored then ^ = 
except when j = (which corresponds to a Jacobian in- 
dependent of /i) . Hence the determinant of the Jacobian 
is a constant. This in turn implies that chiral restoration 
in the average sense means that spatially averaging must 
be equivalent to a uniform averaging over the internal 
space, precisely as one would expect intuitively. 

To proceed further, let us consider a broader class 
of order parameter: G (i7(^(x))) S'(x), where G is a 
chiral order parameter of the class considered previ- 
ously, namely a function of U only and 5' is a chi- 
ral scalar of U (but may be non-local e.g., ^(x) = 
C/^(x)[/(x -f- d), with d a vector- valued fixed point, or 



may potentially also depend on derivatives or integrals 
e.g., / (i'^?/Tr[[/'^(x)?7(x+y)] exp(— ^). It is easy to see 
that a necessary condition for chiral restoration in the 

average sense is 
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G(Ai(x))5(x) dV 



(11. 



This follows since if the system is chirally restored in the 
average sense, chiral rotations do not affect the integral 
on the left-hand side. Thus, we can chirally average with- 
out affecting the result. On the other hand, S'(x) is unaf- 
fected by chirally averaging. Moreover, when averaging 
over chiral space for G, it does not matter where in chiral 
space one starts, one gets the same result; thus the chiral 
average is independent of the value of ?7(x) and hence 
of x. One can therefore remove the chirally-averaged G 
from the spatial average yielding Eq. (|IL8)) . 

Since G is an order parameter, the first factor on the 
right-hand side is necessarily zero and the integral on the 
left-hand side must vanish. Thus 
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G (C/(Ai(x))) S'(x) dV 



G[^l)S{^{^l))J{^i)d^l = Q, 



(11.9) 



where the second form follows from change variables of 
integration to /x; as before we can always cast the integral 
into this form even if the mapping is not one-to-one with 
one region in /x corresponding to more than one region 
in X by summing over the contributions of the Jacobian 
time S for the different regions in x. At this stage we 
can simply repeat the argument used previously for F to 
show that chiral restoration in the average sense requires 
that S{x.{^))J{ii) must be a constant independent of /i. 
On the other hand, we have already shown that J is 
a constant. Thus, we conclude that S also must be a 
constant. 

Note that this is a very strong constraint. Since chiral 
restoration means that all chiral order parameters must 
vanish and since the previous argument holds for any 
chiral scalar operator, we conclude that in order for chiral 
restoration in the average sense to occur, all chiral scalars 
must be uniform in space. This completes the first part 
of the proof. 



B. Spatial uniformity for all chiral singlets is 
incompatible with chiral restoration 

The next step is to show that in Skyrme-type mod- 
els it is not possible for all chiral singlets to be spatially 
uniform while simultaneously having chiral symmetry re- 
stored (in the average sense). We do so via an indirect 
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proof: we will assume that the system is in a chirally re- 
stored phase and then show that if all chiral singlets are 
spatially uniform, there is a mathematical inconsistency. 

A key point about chirally restored phases that was 
stated in the previous section is that spatially averaging is 
equivalent to averaging uniformly over the internal space. 
This in turn means that every value of U occurs at some 
point in space (and indeed in a crystal at least once per 
unit cell). Thus, there must exist four points xq, xi, X2 
and X3 with the property that 



(7(xo) = / 

f/(Xl) = ITl 
{/(X2) = iT2 
?7(X3) = ITS 



(11.10) 



where / is the two-dimensional identity matrix. 

To proceed further, we introduce three matrix-valued 
functions: 



0i(x) = C/(x)C/+(x + xi-Xo) 

02(x) = ?7(x)t/t(x + X2-Xo) 

03(x) = C/(x)C/t(x + X3-xo) 



(11.11) 



Next observe that under chiral transformations 0j(x) 
transforms according to 6'j(x) — )• L6'j(x)L^ where L is 
the SU(2) matrix generating left chiral rotations. Note 
that by construction 



;(xo) 



(11.12) 



The analysis in this section is based on the chiral prop- 
erties of quantities constructed from 9j. 

It is obvious from the transformation properties that 
Tr {9j{x.)) is a chiral invariant. Since we are in a chirally 
restored phase in the average sense (by hypothesis), this 
requires Tr(6'j(x)) to be translationally invariant. How- 
ever, by construction Tr(0j (xo)) ~ 0, from which one 
sees that Tr(6'j(x)) = for all x. Since the 6j{x.) are 
generically traceless SU(2) matrices, they can be written 
in the following form: 



,(x) 



n'^{x)Ta where n°(x) 



.Tr(6'j(x)T") 
' 2 ■ 
(11.13) 

For our purposes, it is useful to consider the n" as being 
components of three distinct unit isovectors each labeled 
by j. Thus nj(x) = (n](x), n^(x), n^(x)) so that 6j{x.) = 
—inj{x) ■ T where arrows will indicate vectors in isospace 
and hats unit vectors in isospace. The fact that fij(x) 
has unit norm follows from the fact that 9j{x.) € SU(2). 

A very useful set of chirally invariant quantities can be 
constructed from the hj{x.): 



iy,;^(x) = ri,(x)-ri,(x + d) 

Tr((?,(x)0j(x + d)) 



where d is a vector-valued parameter. It is quite straight- 
forward to show that for any values of i,j and d, 
is, in fact, a chiral singlet. Note that given the result 
demonstrated in the previous section, if the system is in 
a chirally-restored phase, then W^{x) must be indepen- 
dent of X for fixed ij and d. 

Clearly, by choosing different values of d, one can 
probe spatial correlations among the Li. One important 
case is where d = 0. For notational convenience we will 
introduce a special symbol for this: 

o«(x,.<,x,.IifflMMi±i!)). ,„,5) 

It is trivial to see from Eq. (jll.lSP and Eq. (|II.12p that 



;(xo) = S,j 



(11.16) 



By hypothesis, the system is in a chirally-restored phase 
and thus Qij (x) is a chiral singlet independent of x. Since 
Qij is simply the inner product of a set of three unit 
isovectors, the most general way for it to be independent 
of position is if the set of h at one point are related to 
the set at another by the most general inner product 
conserving transformation, an element of S0(3). Thus 
we conclude that 



hj{x} = y(x)?ij(xo) with^(xo) e S0(3) 



(11.17) 



where R (x) is the same for all j. Note that assuming 
U describes a configuration, which is chirally restored in 
the average sense, ^(x) is completely determined by the 
chiral field U (x) . Recall that 9j (x) transforms accord- 
ing to 0j{x) — !■ L6j{x)L^ and are thus in the irreducible 
(^,^) chiral representation. As such they are chiral or- 
der parameters and by definition their spatial averages 
must vanish in a chirally-restored regime. This in turn 
means the spatial average of the Uj must be zero. Using 
Eq. ([lL2l) this means that in a chirally restored phase 

/ dl/1?(x)n,(xo) =0 (11.18) 

for all J. However since the fij form a complete and 
orthonormal basis, this is possible only if the integral of 
the matrix itself vanishes: 

f dV^{x) = (11.19) 

The hypothesis that the system is in a chirally-restored 
phase (in the spatially-averaged sense) requires that all 
chiral singlet observables are independent of position. 
Since Wj^(x) is chiral singlet, it must be independent 
of X. Using Eq. pi.l7P and the definition of W from 
Eq. (jII.14[) we see that 

Wt^{x) = fif (xo)'i?^(x)1^(x + d)n,(xo) (H.20) 



(II 14) Note that W7^(x) can be expressed as the matrix 
R "^(x) i? (x + d) evaluated between the vectors 7ij(xo). 
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However, since these vectors form a complete basis and 

are independent of x, the only way for WfAx) to be in- 

^ y ^ y 

dependent of x is if the matrix R ^(x) _R (x + d) itself is 
also independent of x: 

^'^(x)^(x + d) = W(d + xo) . (11.21) 
But since the transpose of a rotation matrix is its inverse, 

^(x + d) = ^(x)^(d + xo) . (11.22) 

Note that x and d are arbitrary three-dimensional vec- 
tors in space. Thus the form of Eq. (jll.22p holds under 
replacements of x and d by other vectors. Consider in 
particular the replacements x — ^ d -I- Xq and d — )- x — Xq. 
Thus i? (x + d) = i? (d -I- xo) i? (x). Equating the two 
ways of writing i? (x + d) and relabeling d + xq as y 
yields 



[^(y), 



— ^ , 

i?(x)] 







(11.23) 



Note that since d was arbitrary so is y and Eq. (jll.23p 
holds for all values of x and y. Eq. (|II.23p strongly 
constrains the nature of R (x) : R at any points com- 
mutes with R at any other point while general three- 
dimensional rotations do not. Indeed, the only way that 
Eq. pi.23p can be satisfied for all points is if the rota- 
tions at every point in space are all about a common 
spatially independent rotation axis in isospace. Let us 
denote the unit vector in the direction of this axis n^ot ■ 
It is a general property of rotation matrices in three di- 
mensions that the axis of rotation is an eigenvector with 
an eigenvalue of unity. Since firot is independent of x it 
is also an eigenvector of the spatially-averaged rotation 
matrix again with an eignenvalue of unity: 



n-rot 



(11.24) 



However, Eq. (jll.24p is clearly inconsistent with 
Eq. ((TLT9)) : the spatially -averaged matrix cannot simul- 
taneously be zero and have an eigenvalue of one. Note 
that both Eq. pi.24p and Eq. pi.l9p are direct mathe- 
matical consequences of our assumption that the system 
is in a chirally-restored phase (in the spatially-averaged 
sense). Equation pi.lQp followed quite directly from this 
as by definition all chiral order parameters, including 
9j, must vanish upon spatially averaging in this phase. 
Equation pi.24[) followed from the requirement proved in 
the previous subsection that all chiral singlets, including 
W^, must be spatially uniform in such a phase. Since 
these two are inconsistent, we have shown that the as- 
sumption underlying them — that the system is in a phase 
which is chirally restored in a spatially averaged sense — 
cannot be correct. Moreover, since the analysis was done 
for an arbitrary field configuration for the entire class of 
Skyrme-type models we have proved that chiral restora- 
tion in the average sense cannot occur in models of this 
class. 



Ultimately, this result should not be surprising. The 
only natural way to conceive of a situation in which chi- 
ral order parameters to be non-zero locally but zero un- 
der spatial averaging is for spatial averaging to uniformly 
cover all directions in the internal chiral space. This will 
occur if the field configurations map from real space to 
chiral space in such a way that uniform coverage in one is 
mapped to uniform coverage in the other. However, the 
internal chiral space is curved while the physical space is 
fiat and hence such a mapping is not possible. 



III. A NO-GO THEOREM FOR LARGE QCD 

In this section, we show that the result derived in the 
previous section for the Skyrme model holds for large Nc 
QCD itself. The precise statement is that if large Nc 
QCD is in a phase in which the chiral condensate (qq) is 
generally non-zero but varies from point to point, then 
it is not possible for all chiral order parameters to van- 
ish under spatial averaging. The strategy for doing this 
exploits the fact that the proof in the previous section 
holds for any field configurations describable as a nonlin- 
ear signia model. Thus our theorem will be established in 
general provided we can show rigorously that in a puta- 
tive chirally-restored phase with a spatially varying but 
nonzero chiral condensate the expectation values for a 
key set of operators are reducible to those of a nonlinear 
sigma model. 

To do this we focus on the scalar-isoscalar and 
pseudoscalar-isovector quark bilinears qq and q~^^^q. 
These operators transform into one another as members 
of a (^,^) representation of the SUl(2) x SUr(2) chi- 
ral group. Suppose that the system is in some known 
quantum mechanical state. Let us combine the expecta- 
tion values of these operators in this state into a single 
two-dimensional matrix- valued function: 

... , (g(x)g(x)) / + (g(x)^7^g(x) ■ ^) .^^^.^ 

v{x) = m.i 

where / is the identity matrix and the normalization fac- 
tor, {qq)viLc is included so that in the vacuum state V 
simply becomes the identity matrix. Note that under chi- 
ral transformations on the operators making up V means 
it transforms in precisely the same way as U does in the 
Skyrme model. One can construct order parameters from 
V. Since we are at large 



{F{V{^)))^F{{V{^))) , 



(in.2) 



where the angle brackets indicate quantum expectation 
values in the state. Thus we can treat the issue of chiral 
restoration in the spatially-averaged sense as if V were a 
classical function. 

To proceed further let us introduce a chiral scalar func- 
tion given by the norm of V: 



u(x) 



Tr(Ft(x)F(x)) 



(in.3) 
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By construction, t;(x) is real and nonnegative. For all 
points in space where w(x) ^ one can define 

"W-^ (,,,,4) 

Note that this construction ensures that [/(x) e SU(2). 
Thus provided that w(x) is non-zero throughout space 
and the expectation values vary smoothly through space, 
we have an associated smoothly varying SU(2) matrix 
[/(x), precisely as in nonlinear sigma models such as 
the Skyrnie model. Note that the argument in subsect. 
Ill Al does not depend on the dynamics but merely on 
the transformational properties plus assumptions about 
smoothness. Thus we conclude that i/ the system is in 
a phase which is chirally restored in the average sense 
while simultaneously having w(x) ^ everywhere, all 
chiral scalars must be spatially uniform. Note moreover 
that u(x) is a chiral scalar itself and therefore must itself 
be a constant in any putative chirally restored phase with 
V everywhere non-zero. 

Next, we need to eliminate the possibility in large 
QCD of v(yi) vanishing at some points in space while 
being in a phase which has a spatially varying chiral 
condensate while being chirally restored in the spatially- 
averaged sense. To proceed let us consider what happens 
if V is distorted slightly away from its physical value: 

V^(x) = V^(x) + A AV^(x) (III.5) 

where A is a small parameter and AX^ is an arbitrary 
smooth function with the property that it is non-zero in 
the neighbor of any point where V vanishes. We will use a 
tilde to denote a quantity perturbed away from its physi- 
cal value due to the use of V in place of V . Note that the 
change in V due to this distortion is small and smooth. 
So the changes in spatially-averaged chiral order param- 
eters obtained from V are expandable as a Taylor series 
in A. Thus j/ the system was in a putative phase, which 
was chirally restored in the average sense, then spatially- 
averaged chiral order parameters computed with V would 
be of order A ; V would describe a nearly chirally restored 
regime. 

To proceed, we exploit the fact that V is constructed to 
be nonzero everywhere, implying that \J is well-defined 
everywhere. We will describe any putative nearly chi- 
rally restored regime associated with V using \J . Since 
the supposed regime is nearly chirally restored, the ar- 
gument in subsect. Ill Al goes through up to corrections 
associated with the fact that it is only nearly restored. 
Thus we see that in such a regime all chiral singlets must 
be nearly constant in space: s(x) = sq + A(5s(x), where 
So is a constant and A(5.5(x) describes the fluctutations 
away from it. Since w(x) is a chiral scalar we conclude 
that u(x) — vq + A(5z;(x). On the other hand, by con- 
struction at any point in space, 'D(x) can be expressed as 
a Taylor expansion in A : w(x) = 'y(x) -I- A Aw(x). Equat- 
ing these two forms yields 

v(yi) = Wo + A ((5w(x) - Aw(x)) . (III.6) 



This implies that as A -^0, u(x) approaches a constant 
value. However, a constant value of w(x) is inconsistent 
with V going to zero at some points while generally be- 
ing non-zero as a result of the non-zero chiral conden- 
sate. Thus z/ a phase which is chirally restored in the 
spatially averaged sense exists and has a spatially vary- 
ing chiral condensate, it must have a constant non-zero 
w(x). Therefore, it is effectively reduced to a nonlinear 
sigma model and the proof in the preceding section ap- 
plies. This completes the demonstration. 



IV. SKYRMIONS ON A HYPERSPHERE 

We have shown that chiral symmetry restoration in 
the spatially-averaged sense is not possible for either the 
Skyrme model or for large QCD. However, two im- 
portant issues arise regarding the Skyrme model. As was 
noted in the introduction observed chiral restoration in 
the spatially-averaged sense is observed for the Skyrme 
model on the hypersphere d, llGl - flsl [2l] | . 

The first issue is quite straightforward. One might 
worry that since chiral symmetry does get restored in 
the average sense for Skyrmions on the hypersphere this 
might indicate there is a flaw in the proof that Skyrmions 
cannot have chiral restoration at finite density. Of course, 
it need not indicate a flaw since the geometry for which 
situation for which the no-go theorem was derived — 
extended nuclear matter in flat space — is different from 
the curved space of the hypersphere. Nevertheless, it is 
important to pin down how the proof can be valid for 
three-dimensional flat space but fail for the hypersphere. 

The second issue concerns intuition. As noted in the 
introduction, the approximation of a Skyrmion crystal 
by a single Skyrmion on a hypersphere is ad hoc. The 
notion, though, was that the main qualitative effect of 
putting a Skyrmion into a crystal is to restrict the space 
over which it can spread. This effect is certainly also 
present with a single Skyrmion on a hypersphere. To 
the extent that qualitative questions such as whether or 
not chiral restoration (in a spatially-averaged sense) oc- 
curs is a generic feature which depends on whether a 
Skyrmion is confined to a sufficiently small volume, one 
might think that a Skyrmion on a hypersphere would 
be a good way to discover this. However, this is clearly 
wrong. Chiral restoration in the average sense occurs on 
the hypersphere but not for crystals in flat space. It is 
important to understand why the intuition gained from 
the hypersphere fails. 

Before addressing these issues, let us briefly review 
the geometry of the hypersphere. This can be charac- 
terized as a three-dimensional curved surface in a four- 
dimensional Euclidean space with fixed radius: 

+ + = (IV.l) 

where L is the radius of the hypersphere. It is useful 
to parameterize this in terms of three angular variables. 
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< ^, 6* < TT and < (/) < 27r with 

X = L sin /.( sin 9 cos ( 
y = L sin ^ sin 9 sin q 
z = L sin /i cos 9 
w = L cos /J, 



(IV.2) 



We denote the metric tensor for this geometry as g^^. 
With parameterization used here it is diagonal and the 
diagonal matrix elements are: 



5mm = 



r 2 • 2 

L sm fl 



(IV.3) 



(700 = sin^ /i sin^ . 

The volume element is given by 

dV = y/dctgd^i d9 dcj) = sin^ ^l dfj. sin 9 d9 d(t> (IV.4) 

Consider the configuration 

[/o(m,^,</') = 

cos(^) + sin(//) cos(6') sin(/x) sin(6') exp(— i0) 
sin(/x) sin(6') exp(i0) cos(/x) — sin(/Lt) cos(6') 

(IV.5) 

It is straightforward to verify that this configuration has 
winding number unity and that any chiral order param- 
eter constructed from Uq vanishes when integrated uni- 
formly over the hypersphere. Thus the configuration 
corresponds to chiral restoration in the spatially aver- 
aged sense on the hypersphere. The explicit example of 
Uq shows that the hyperspherical geometry does allow 
configurations with chiral restoration in the spatially- 
averaged sense. The theorem proved in Sect. HIl showed 
such a configuration — at least for the case of flat space. 
This means either the proof of the theorem is wrong or 
some aspect of the proof holds for flat space but fails for 
the case of the hypersphere. 

Fortunately, it is easy to see how the proof of the 
no-go theorem can hold for flat space and not for the 
hypersphere. At several critical point points in the 
derivation the fact that the space is flat plays a critical 
role. Consider as an example the relation i? (x + d) = 
R (x) R (d + xo) from Eq. ((IL22|; th is relat ion is central 
in the analysis leading to Eq. pi.23p which is at the core 
of the theorem. However, as written this result is not 
meaningful on the hypersphere: there is no notion of lin- 
early adding the vectors associated with two points in the 
space to obtain the vector associated with a third point 
in the space. Thus the notion of the point x + d is simply 
ill-posed on the hypersphere and the proof developed in 
flat space does not go through. 

Let us now turn to the issue of the failure of the intu- 
ition that a single Skyrmion in the hypersphere should be 
qualitatively similar to the Skyrme crystal as the princi- 
pal effects of both should be to confine a Skyrmion to a 



limited spatial region. The intuition can break down for 
one of two reasons: either limiting the space in which a 
Skyrmion extends by a single Skyrmion in a compact ge- 
ometry is not qualitatively similar to limiting it by plac- 
ing in a crystal in fiat space, or because the hypersphere 
is an atypical compact geometry. It is instructive to un- 
derstand which is the cause. There are good reasons to 
suspect that it is due the atypical properties of the hy- 
persphere. The hypersphere has a much higher symmetry 
than typical geometries one can consider. It is plausible 
that these symmetries rather than generic properties are 
responsible for the chiral restoration (in the spatially- 
averaged sense) seen for Skyrmions on the hypersphere. 
This becomes particularly plausible when one considers 
the internal space associated with the field U . It can be 
parameterized as 



U ~ si + ip- f 



(IV.6) 



where s and p can be extracted from U using the relations 



Tr(C/) 



Pj = 



(IV.7) 



Note that by construction s and p satisfy the constraint 



pI 



pI = ^ 



(IV.8) 



This constraint in the internal geometry in Eq. (jlV.Sp is 
of the same form as the constraint in Eq. PV.1[) describ- 
ing the hypersphere. Indeed it is easy to see that the 
configuration in Eq. pV.5[) which restores chiral symme- 
try on the hypersphere (in the spatially averaged sense) 
has the property that it maps points on the hypersphere 
to the analogous points in the SU(2) matrix U: 



(IV.9) 



Thus averaging over the hypersphere uniformly automat- 
ically leads to averaging uniformly over the internal space 
forcing all chiral order parameters constructed from U to 
vanish. 

One way to test whether the observed chiral restoration 
seen for a Skyrmion on a sufficiently small hypersphere 
is due to the special geometrical properties of the hy- 
persphere is to ask what happens to spatially-averaged 
chiral order parameters in this regime if the geometry 
is distorted slightly away from a hypersphere. If chiral 
symmetry is generically restored for single Skyrmions in 
small compact spaces, these spatially averaged order pa- 
rameters should remain zero. If instead chiral restoration 
is a result of special features associated with the hyper- 
spherical geometry one would expect these to become 
non-zero. 

To make this analysis concrete, we need to impose dy- 
namics associated with a particular variant of the Skyrme 
model. Here, we will use the simplest one which contains 
field U and up to four-derivative terms: 

f2 2 

£ = ^^Tr[d^,U''d''U] + jTr[U''df,U, U"^ d^Uf , (IV.IO) 
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where f^^ ~ 93MeV is the pion decay constant and e is 
a dimensionless parameter necessary to stabihze the soh- 
ton. In the remainder of this paper, we will use dimen- 
sionless units by setting /^r and e to unity. One can use 
dimensional analysis to reinsert factors of and e at the 
end of the problem, if desired. To do so, simply multiply 
lengths by 2\/2e//7r and energies by V2eU uM- For the 
purposes of determining whether or not chiral symmetry 
is restored, these rescaling factors are irrelevant. 

For the case of an unperturbed hyperspherical geom- 
etry, the minimum energy winding number unity con- 
figuration is a hedgehog localized in one region of the 
hypersphere, providing the radius L is smaller than the 
critcal value of [l3|. For L < V2, the minimum en- 
ergy configuration is given by Uq in Eq. PV.SP or config- 
urations obtained from it by a global SU(2) rotation in 
internal space. Our strategy is to start with the lowest 
energy configuration for a value of L well into the restored 
phase (for concreteness we take L = 1). We then con- 
sider a small perturbation on the geometry, compute the 
shift in the configuration due to this to lowest order in 
the perturbation, and use this perturbed configuration 
to compute the spatially-averaged chiral order parame- 
ters at first order in the shift. If chiral restoration in the 
spatially-averaged sense were a generic feature seen when 
a Skyrmion is confined by a compact geometry of suffi- 
ciently small size one expects spatially averaged order 
parameters to remain zero. 

The most general perturbation of the hypersphere can 
be formulated as a shift in the contravariant metric per- 
turbation away from that of the hypersphere parameter- 
ized in terms of an overall expansion parameter A : 

9'Ut-9l + ^l'' ^ (IV.ll) 

where 7*'' is dimensionless. We will consider a very simple 
class of perturbations to the geometry: those in which 
all matrix elements of 7 vanish except the {i = 1, j = 1) 
component. Further, we will take 7^^ to be a sinusoidal 
function of ^ onl,y which ensures that the perturbation 
does not induce discontinuities in derivatives: 

f J = ^Ji sin(p/i) , (IV.12) 

where p is a non-zero integer. Of course, this is a very 
restricted class of perturbation. If spatially-averaged chi- 
ral order parameters remain zero with it, it would tell us 
very little since this might be a result of the highly con- 
strained and symmetric form of the metric perturbation. 
However, if they do not remain zero, then it is sufficient 
to show that the hyperspherical geometry does not be- 
have generically. 

We now solve the Skyrme model for the lowest energy 
unit winding number solution with this geometry working 
to first order in A. To do so we consider configurations 
of the most symmetric sort consistent with a winding 
number of unity. These are hedgehog configurations of 



the form 

(sin(6') cos((/)) \ 
sin(6') sin(0) 
cos{e) J 

with /(O) = and /(tt) = n . 

(IV. 13) 

We note that for the unperturbed geometry, Uq, the con- 
figuration in Eq. (jIV.5[) associated with chiral restora- 
tion is of this form with /(/i) — ji. For the perturbed 
case, we exploit the fact that the most symmetric form 
remains consistent with the topological constraints. It 
is guaranteed that the configuration of this form, which 
minimizes the energy, will also be a solution of the full 
Euler-Lagrange equations. As a general rule, a solution 
of the Euler-Lagrange equations obtained from a highly 
symmetric ansatz need not correspond to a global min- 
imum of the energy. However in this case it should, as 
in the unperturbed case where the hedgehog is known to 
be a global minimum for that problem and this problem 
differs from that only to first order in A. Next we take 
the ansatz in Eq. (jIV.13[) using 

f{^i) = ^l + \5f{^L) with <5/(0) = <5/(^) = . (IV.14) 

Using this form, we compute the energy to second order 
in A (noting that first-order perturbations in the field 
correspond to second-order perturbations in the energy) 
and vary this form to get differential equations for 5f 
which are then solved numerically. Using these numeri- 
cal solutions the spatial averages of various chiral order 
parameters are computed. 

A useful set of chiral order parameters to consider are 
local functions constructed from U or, equivalently. from 
s and p from Eq. pV.7p . We will focus here on three 
representative samples from this list: 

^_d±4±ii and (IV.15) 
4 4 

It is clear that s and are chiral order parameters since 
s — >■ — s is a chiral transformation. Similarly, it should 
be clear that 2| P1+P2+P3. ^ chiral order parame- 
ter since in a chirally restored phase one necessarily has 
(s^) = {p\) = {p\) = and thus the quantity vanishes. 
In Table HI we show the spatial averages of these order 
parameters computed to leading order in A. 

It is noteworthy that not all entries in this table are 
zero. Thus making small perturbations away from the 
hyperspherical geometry yields non-zero order parame- 
ters. From this we see that chiral restoration seen for 
Skyrmions on the hypersphere is not generic but rather is 
a consequence of the special geometric properties of the 
hypersphere. Of course, this geometry has no physical 
significance and was used for ease of computation. While 
it was hoped that a calculation in this simple situation 
would give useful physical intuition about the more com- 
plicated situation of Skyrmion crystals, it appears that 
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Table I: Spatially-averaged chiral parameters at order A. 
These are computed using geometries of the perturbed hyper- 
sphere given in Eq. (|IV.1H) . The configurations used to com- 
pute these are minimum energy configurations of the hedge- 
hog form based on the Lagrangian in Eq. I|IV.10[) . 



metric perturbation 


s 


3s^ P1+P2+P3 
4 4 




7^^' = S'^ sin(/i) 





-0.0954931 A 





^ §]i sin(2/i) 








-0.0314619 A 







-0.0197049 A 





7*^' = 5'^ 5'^ sin(4/i) 








-0.00699138 A 


= 5^1 sin(5/i) 





0.0156628 A 






the opposite is true with regard to chiral symmetry. In- 
deed, the special properties of the geometry which makes 
analytic calculations of the small radius (high density) 
phase simple also make the intuition totally unreliable 
even for qualitative issues associated with chiral symme- 
try breaking and its possible restoration in the average 
sense. 



V. DISCUSSION 

The question of whether dense, cold nuclear matter 
is chirally restored above some critical baryon chemical 
potential remains a problem of central importance to nu- 
clear physics. It remains unanswered. The answer is also 
potentially unknown in the simpler case of large Nc QCD. 

This paper dealt with a complication at large N^- the 
fact that nuclear matter is likely to crystallize and thus 
the breaking of chiral symmetry and of translational sym- 
metry may get entangled. One could imagine a situation 
in which the breaking of translational symmetry can in- 
troduce non-zero but spatially-averaged chiral order pa- 



rameters which integrate to zero over space. To deal with 
this complication, it was suggested here that the natural 
focus should be on spatially-averaged chiral order param- 
eters which by construction are insensitive to details of 
how translational symmetry is broken. We demonstrated 
here that while it is possible for some spatially-averaged 
chiral order parameters to vanish, it is not possible for all 
of them to if the chiral condensate is generally non-zero 
but spatially varying. Thus chiral symmetry restoration 
in a spatially averaged sense is not possible at large Nc 
unless the chiral condensate is zero everywhere. 

This result is of some significance in connection to 
the question of whether chiral restoration occurs at high 
baryon density in large Nc QCD. It has been argued on 
the basis of Skyrme models that chiral restoration does 
occur at sufficiently high density. The fact that at high 
densities the lowest energy configurations in Skyrme crys- 
tals have a vanishing spatially-averaged chiral condensate 
has be taken as support for chiral restoration Q. As 
shown here, however, there must be other chiral order 
parameters which are non-zero upon spatial averaging. 
It has similarly been argued that the vanishing of chiral 
order parameters upon spatial averaging for a Skyrmion 
on a sufficiently small hypersphere has also been taken as 
evidence for chiral restoration at high density Q. How- 
ever, as shown here this was an artifact of the hyperspher- 
ical geometry, and is thus unconnected to the question of 
what happens in crystals in flat space. 
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